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Abstract 

The mesonic excitations and s-wave tt-tt scattering lengths at finite temperature are 
studied in the two-flavor Polyakov-Nambu-Jona-Lasinio (PNJL) model. The masses 
of tt meson and a meson, pion-decay constant, the pion-quark coupling strength, and 
the scattering lengths ao and at finite temperature are calculated in the PNJL model 
with two forms of Polyakov-loop effective potential. The obtained results are almost 
independent of the choice of the effective potentials. The calculated results in the PNJL 
model are also compared with those in the conventional Nambu-Jona-Lasinio model 
and indicate that the effect of color confinement screens the effect of temperature below 
the critical one in the PNJL model. Furthermore, the Goldberger-Treiman relation and 
the Gell-Mann-Oakes-Renner relation are extended to the case at finite temperature 
in the PNJL model. 
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1 Introduction 



QCD thermodynamics and phase diagram, especially about the restoration of the chiral 
symmetry and the deconfinement phase transition which are expected to occur in ultra- 
relativistic heavy-ion collisions [1-8] or in the interior of neutron stars [9-12], has been a 
subject of intense investigation in recent years. One significant aspect to investigate the 
restoration of the chiral or axial symmetry and the deconfinement phase transition is to 
study the variation of properties of particles propagating in hot and/or dense medium [13- 
15]. In this work, we focus on the influence of a hot medium on the properties of light 
pseudoscalar (n) and scalar (cr) mesons, and n-n scattering lengths. Special attentions are 
paid to their dramatic variations near the regime where the chiral phase transition and the 
deconfinement phase transition occur. We expect to extract the signals of phase transition 
from mesonic excitations and ir-ir interactions in the hot medium. 

A promising phenomenological approach to study the low-energy processes involving the 
pseudoscalar and scalar mesons at zero temperature and finite temperature is the Nambu- 
Jona-Lasinio (NJL) model [16-21]. The most important advantage of the NJL model is 
that it introduces a mechanism of the dynamical breaking of chiral symmetry (due to the 
quark-antiquark condensate). However, the NJL model has its disadvantage, which is the 
lack of the description of color confinement. To include some effects of color confinement, 
a Polyakov-loop improved Nambu-Jona-Lasinio (PNJL) model has been developed recent 
years [22-30]. In the PNJL model, the Polyakov-loop as a classical field couples to quarks 
and thus suppresses the contributions from wrong degrees of freedom (color non-singlet) 
to the thermodynamics below the critical temperature. Therefore, the introduction of the 
Polyakov-loop represents some aspects of the color confinement, at least on the level of 
statistics [29] . The validity of the PNJL model has been confirmed in a series of works by 
confronting the PNJL results with the lattice QCD data [27, 28, 31-33]. The phase structure 
and thermodynamics in the PNJL model have recently been explored extensively [29, 30, 34- 
46] , and the impact of Polyakov-loop dynamics on the chiral susceptibility or quark number 
susceptibility [34,35], QCD critical endpoint [40,41] and critical surface [42], and the color 
superconductivity phase transition [28,47-49] have attracted lots of interests. Furthermore, 
fluctuations beyond the mean field approximation have been included in the PNJL model [50, 
51], and the PNJL model has also been extended to the regime of imaginary chemical 
potential [52-55] and + 1 dimensions [56], and applied to analyze the flavors of quark- 
gluon-plasma [57] and the isentropic trajectories on QCD phase diagram [58]. 

The properties of pseudoscalar and scalar mesons at finite temperature for two [14] and 
three [15] flavor systems have also been investigated in the PNJL model. In Ref. [14], The 
mesonic correlators and spectral functions for ir and a mesons were obtained. It was found 
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that the n-a degeneracy in the chiral symmetry restored phase was still satisfied after cou- 
pling quarks to the Polyakov-loop and the role of ir meson as Goldstone boson was also 
confirmed in the PNJL model. It was also found that, although the PNJL model can not 
cure the problem of the conventional NJL model as for the unphysical width of the a me- 
son, the PNJL results on the decay width improved slightly the NJL ones [14]. In order 
to further study the broken chiral symmetry and its restoration in the mesonic sector in 
the PNJL model which makes the investigation of the interplay between the restoration of 
chiral symmetry and the deconfinement phase transition possible, it is necessary to study the 
Goldberger-Treiman relation [59] and the Gell-Mann-Oakes-Renner relation [60] which are 
direct results due to the chiral symmetry breaking. Furthermore, one of the most fundamen- 
tal hadronic processes of QCD at the mesonic level, the pion-pion scattering, n + n — > n + n, 
at finite temperature, which provides a direct link between the theoretical formalism of chiral 
symmetry and experiment, also deserves to be investigated. In this work, we will then study 
the problems mentioned above in the PNJL model. 

The paper is organized as follows. In Sec. II we simply review the formalism of the two 
flavor PNJL model. In Sec. Ill we discuss the mesonic excitations at finite temperature in 
the PNJL model. The dependence of the pion-decay constant, pion-quark coupling strength, 
and the relation between the mass of a meson and that of tt meson on the temperature are 
studied in the PNJL model. We also extend the Goldberger-Treiman relation and Gell- 
Mann-Oakes-Renner relation to a formalism which is appropriate at finite temperature. In 
Sec. IV we study the s-wave tt-tt scattering lengths in the PNJL model and compare the 
results in the PNJL model with those in the conventional NJL. Finally, in Sec. V, we give 
a summary and conclusions. 



2 The PNJL model 

The Lagrangian density for the two-flavor PNJL model is given as [27] 

Cpnjl = 4> (hn D>l ~ ™o) 4> + G (inp) 2 + (#7 5 r^) 2 j 

-W($[A],<5L4],T), (1) 

where ip = (ip u ,i J d) T is the quark field, 

D» = d t *- %A* with A<* = 5% A , A = ^4°y = -iM- (2) 

The gauge coupling g is combined with the SU(3) gauge field A%(x) to define A^(x) for 
convenience and A a are the Gell-Mann matrices in color space. rh = diag(m u , m^) is the 
current quark mass matrix. Throughout this work, we take m u = = m , assuming the 
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isospin symmetry is reserved on the Lagrangian level. The four-fermion interaction with an 
effective coupling strength G for scalar and pseudoscalar channels has SU V (2) x SUa(2) x 
Uv(l) symmetry, which is broken to SUv(2) x Uy(l) when m 7^ 0. Here r a (a = 1,2,3) in 
the Lagrangian density (Eq. (1)) are Pauli matrices in flavor space. 

The U ($, l», T) in the Lagrangian density is the Polyakov-loop effective potential, which 
controls the Polyakov-loop dynamics and can be expressed in terms of the trace of the 
Polyakov-loop $ = (Tr c L)/iV c and its conjugate $ = (Tr c V)/N c . Here the Polyakov-loop L 
is a matrix in color space, which can be explicitly given as [27] 

rP 



L (x) = V exp 



dr A4 (x, t) 



exp [if3A 4 ] , (3) 



where (3 = 1/T is the inverse of the temperature. The Polyakov-loop effective potential has 
the Z{3) center symmetry like the pure-gauge QCD Lagrangian. When the temperature is 
lower than a critical value (T ~ 270 MeV in pure gauge QCD [27]), the value of $ (and 
$) which minimizes the Polyakov-loop effective potential is zero, meaning that the phase is 
color confined and has the Z{3) symmetry. However, when the temperature is above the 
critical temperature T , $ develops a nonzero value which minimizes the effective potential 
and the system is transited from a Z{3) symmetric, confined phase to a Z{3) symmetry 
broken, deconfined phase. The temperature dependent Polyakov-loop effective potential is 
chosen to reproduce the lattice data for both the expectation value of the Polyakov-loop [61] 
and some thermodynamic quantities [62]. In the PNJL Lagrangian in Eq. (1), the coupling 
between the Polyakov-loop and quarks is uniquely determined by the covariant derivative 

In previous works, two possible forms for the Polyakov-loop effective potential have been 
well developed. Following our previous work [29], we denote them as W po i($, $,T) and 
Wimpt^, T), respectively. The former is a polynomial in $ and $ [27] and the latter is 
an improved effective potential in which the higher order polynomial terms in $ and $ are 
replaced by a logarithm [28]. Both the effective potentials are taken in our work to investigate 
whether our results depend on the details of the Polyakov-loop effective potential. These 
two effective potentials have the following forms 

with 

b 2 (T) =a + ai (^j + a 2 (^j 2 + a 3 (?f) \ (5) 

U imp (<M,T) = _i A ( T )$ $ + B ( T ) ln |- 1 _ 6 $ $ + 4 ($3 + $ 3)_3($ $ )2] ? (6) 

J. 2 



and 
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with 

A(T)=A + A 1 (fy+A 2 (^y, B(T) = B 3 (^y. (7) 

A precise fit of the parameters in these two effective potentials has been performed to repro- 
duce some pure-gauge lattice QCD data in Refs. [27,28]. The results are listed in Table 1, 
Table 2, respectively. The parameter T is the critical temperature for the deconfinement 
phase transition to take place in the pure-gauge QCD and T is chosen to be 270 MeV 
according to the lattice calculations. 



Table 1: Parameters for the polynomial effective potential U po \ 



a 


ai 


a 2 


a 3 


h 


&4 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



Table 2: Parameters for the improved effective potential IA\ 



A* 


A 1 


A 2 


B 3 


3.51 


-2.47 


15.2 


-1.75 



In the NJL sector of the model, three parameters need to be determined: the three- 
momentum cutoff A, the current quark mass m , and the coupling strength G. In our work 
we employ the zero-temperature values of the quark condensate, pion decay constant and 
the mass of pion to fix the parameters. The obtained results are given in Table 3. 



Table 3: Parameters of the NJL sector of the model and the physical quantities being fitted 



A (MeV) 


G (GeV- 2 ) 


m (MeV) 


\$M V3(MeV) 


U (MeV) 


77V (MeV) 


659.28 


4.773 


5.32 


250.0 


92.4 


139.3 



3 Mesonic excitations at finite temperature 

Before we study the properties of mesonic excitations at finite temperature in the PNJL 
model, the gap equation whose solution provides the constituent mass of the quark should 
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be given. As presented in Ref. [14], such gap equation in the Hartree approximation reads 

+°° r A3. 



m 



= m + 2GTTr V / — ~ — T , (8) 



where the imaginary time formalism is used and the temporal component of the four- 
momentum is discretized, i.e. po = ioj n and u n = (2n+ l)nT is the Matsubara frequency for 
a fermion; m is the constituent mass of the quark; Tr is the trace which operates over Dirac, 
flavor, and color spaces. Here the three-momentum cut-off is employed. After a sum of the 
Matsubara frequencies, Eq. (8) can be written as 

m = m + 2GN f £ / {1 - / [E p - (-iA 4cc )\ - f [E p + (-iA 4cc )]} , (9) 

where E p = (p 2 + m 2 ) 1 / 2 and the summation over the color index can be further written as 

N c 

J2f[E p -(-iA 4cc )\ 

c=l 

N c 

= Y - 

/L^i e /3E pe il3A 4cc _|_ I 
c=l 

= [( e ^p e ^ 4 22 _|_ Ij ( e PEp e iPA 433 _|_ ]_) _|_ ^Ep e i/3A 411 + ^ (^JiE v e if3A 433 + ^ 

+ (eP Ep e il3Mll + l) (eP Ep e i/3A422 + l)] [(eP Ep e i(3Mll + l) (e pEp e ipM22 + 1) {e pEp e ifiAi33 + i)] _1 

= ^1 + ZQe-PE* + 3$e- 2 ^ + e- 3 ^ = (-Ep), (10) 
where the distribution function f£(E p ) in the PNJL model has been given in Ref. [14] 
with another method and we follow their notations. We can find that, when $ = $ = 1, 
f£(E p ) becomes the conventional Fermi-Dirac distribution function. In the same way, the 
summation of the last term in Eq. (9) is 

J2f[E p + (-iA 4cc )} 

c=l 

$ e -/3£ P _|_ 2 $ e -2/3£p + e -3(JE p 

= Nc \ + SQe-PE, + 3|, e -2/3£; p + e -3/3E p = Nc f* ( U ) 
Finally, the gap equation is given by 

m = m + 2GN f N c jf [l - /+(£,) - . (12) 

The gap equation in the PNJL model at finite temperature can also be simply derived from 
the gap equation at zero temperature, which is 

m = m + 8GmN f N c iI 1 , (13) 
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where 



h 



d 4 p 



(2tt)V 



(14) 



To calculate the integral I\ at finite temperature in the PNJL model, we just need to replace 
the integral in p with iT J2 n w J2 C with Po = ~ ^4cc> i- e - 



+00 N c 

iT E ^E 



n=— 00 c= 



d 3 p 



d 3 p 



(2tt) 3 2£ ? 



1 A (2vr) 3 - zA 4cc ) ; 
[1 -/+(^) 



E p 2 



(15) 



We have shown that calculations at finite temperatures in the PNJL model can be sim- 
ply derived from calculations at zero temperature above. Therefore, in the following we 
investigate the mesonic excitations at finite temperature in the PNJL model starting from 
those at zero temperature. We follow the formalism in Ref. [18] and the tt and o mesons 
correspond to the pseudoscalar isovector modes and the scalar isoscalar mode, respectively. 
For the pseudoscalar modes, defining the operators 



T 



V2 



(ri ± ir 2 ) 



(16) 



we can reexpress the four-fermion term in the pseudoscalar channel in the Lagrangian in 
Eq. (1) as 

(ipi^fip) 2 = 2 (ipi^r+ip) (ipi^r-ip) + (^ry 5 T 3 V>) (V^s^) • (17) 

The effective interaction resulting from the exchange of a tt meson can be obtained as an 
infinite sum of loops in the random-phase approximation (RPA) [18] and the leading order 
terms in N c is shown diagrammatically in Fig. 1. 




Figure 1: Schematic representation of the effective interaction for the pseudoscalar modes in 
the RPA, where the double dashed line represents the effective propagator of tt mesons and 
the solid lines are quark lines; the black dots denote the effective coupling between tt meson 
and quarks. Here, Tj = Tj = r 3 for tt°, and Tj = r 1 * 1 , T~ = t t for tt ± . 
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Using the symbols in Ref. [18], the left hand side of the equation in Fig. 1 can be denoted 



as iUij(k 2 ). Summing up all the terms on the right-hand side, we obtain 



(18) 



Comparing Eq. (18) with the equation in Fig. 1, one can find that the mass of it mesons is 
related to the pole of Eq. (18), which is the solution of the following equation [18] 

1 - 2GTI ps (A; 2 ) = 0. (19) 

Furthermore, the coupling strength between n meson and quarks g nqq can be obtained as 

(20) 



2 

9nqq 



\du ps (k 2 y 




ok 2 





Therefore, the information of 7r mesons is included in the pseudoscalar polarization U ps (k 2 ), 
which reads 

-tll ps (k 2 ) = - J ^Tr [t l5 T t tS(k + p)t l5 T,tS(p)} , (21) 

where iS(p) = ij (j> — m) is the propagator of quarks. After calculating the trace in Eq. (21), 
one has 



-ili ps (k 2 ) = AN c Nf J 



d A p 



1 



(2tt)V -m 2 
= 4N c N f I 1 -2N c N f k 2 I(k), 



- 2N c N f k 2 



d A p 



1 



(2tt) 4 {p 2 - m 2 )[{k + p) 2 - m 2 



(22) 



where we have used the function l\ given in Eq. (14) and also defined the function I(k) with 
the same symbols as used in Refs. [18,63,64], i.e. 



4 (p 2 — m 2 ) [{k + p) 2 — m 2 } 



(23) 



Furthermore, we introduce another two functions as done in Refs. [63, 64], which will be used 
in the following: 

d 4 » 1 
p (OA) 

(2vr) 4 {p 2 - m 2 ) 2 [{k + p) 2 - m 2 ]' 1 ' 

(2tt) 4 {p 2 - m 2 ) 2 [{k + p) 2 - m 2 ] 2 ' 1 ' 

Then, substituting the expression of the pseudoscalar polarization in Eq. (22) into Eq. (19), 
we have 

1 - 8GN c N f ih + 4GN c N f k 2 iI(k) = 0. (26) 



K(k) = 
L(k) = J 
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Upon inserting the gap equation (in Eq. (13)) into the above equation, one obtains [18] 

^ + 4GN c Nfk 2 iI(k) = 0, (27) 
m 

whose solution gives the mass of the pseudoscalar mode. The explicit expression for the 
coupling between n meson and quarks can be easily obtained upon substituting Eq. (22) 
into Eq. (20), and, in turn, it reads [63] 

= N c N } I{m^)+I{0)-mlK{m n )- (28) 

In the case of finite temperature in the PNJL model, we need to extend the function I(k) 
in the same way as taken for the function I 1 . When the three momentum is vanishing, i.e. 
k = (uj, 0), I(lu, 0) at finite temperature in the PNJL model is 

°> = "* / Wf E p (J-4E}) I 1 - ^ - WW ' < 29 > 
Then Eq. (27) at finite temperature can be rewritten as 

£ + 4gn < n >< L (S«^) i 1 - it{E - ] - ^ = °- (30) 

Furthermore, we have 

K{u > 0) = *L wr iEW^f t 1 - ~ aw • (31) 

^ °» = 'L S ^3f [i - mEr) - um ■ (32) 

In the same way, properties of a mesons can be extracted from the scalar polarization 
U s (k 2 ), which is 

-iU s (k 2 ) = -J-^Tr[iS(k+p)iS(p)} 

= 4N c N f I 1 -2N c N f (k 2 -4m 2 )I(k). (33) 

Therefore, employing the RPA approximation for the scalar channel in the same procedure 
as for the effective interaction in the pseudoscalar channel, one could determine the mass of 
a meson which is the pole of its corresponding effective propagator, i.e. 1 — 2GU s (m 2 ) = 0, 
explicitly given by 

— + AGN c Nf(m 2 - Am 2 )il{m a ) = 0. (34) 
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The relation between the mass of a meson and that of tc meson could be obtained by 
comparing Eq. (34) and Eq. (27) as 

< = tm ' + <& (35) 

which returns to the relation given in Ref. [18] when the difference between I(m n ) and I{m a ) 
is neglected. 

In the following, we would investigate the pion-decay constant f w at finite temperature 
in the PNJL model with a starting of the definition of f n 

(0\J^(x)\ni) = tkjjv . (36) 

Considering the explicit expression of the left hand side in the PNJL model 



T 



(37) 



we arrive at 

f n = -AiN c g^ qq mI(m^). (38) 

Employing the expression of the pion-quark coupling in Eq. (28) and considering the fact 
Nf — 2 in the present case, one has [63] 

* = + (39) 

and 

/V = 4m 2 - 2 = m 2 r 2 , (40) 

** qq [I(0) + I(m n )-mlK(m n )] 2 1 J 

where we have defined a symbol r as 

= 2/(ttv) (41) 

When the temperature is approaching zero, I(m w ) ~ 1(0), K(m n ) ~ 0, r — > 1. Eq. (40) 
returns then to the quark level version of the Goldberger-Treiman relation at zero temper- 
ature [59]. As we show below, when the temperature is near some critical temperature, r 
deviates from 1 evidently. Furthermore, from Eq. (27), one has the mass of n meson as 

m 2 = _a^o i m 

m 4GN c N f iI(m n ) ' V ; 

One may also combine this equation with Eq. (39). It gives consequently 

2 2 = mom IjjnJ 

W?r ~ G I(0)+I(m w )-mlK(m n )- 1 } 
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And the constituent mass of the quark m is related with the condensate of quark by 

m = —2GNf{uu) + m 

= -2G(i>i>)+m . (44) 

Replacing the constituent mass in Eq. (43) with the quark condensate in Eq. (44), we obtain 

2 

m lfl = -m {ipip)r + — r = -m (^)r |l + ^TD^at ] ■ ( 45 ) 

As mentioned above, at zero temperature r — > 1. Considering the lowest-order contribution 
in mo, one obtains then m%f% — —mofyifj), which is the lowest-order approximation to the 
Gell-Mann-Oakes-Renner relation [60]. 

Before our numerical calculations, we need to give the equations to determine the values of 
the Polyakov-loop $ and its conjugate $. In the mean-field approximation or equivalently the 
Hartree approximation, the thermodynamical potential density for the Lagrangian density 
in Eq. (1) is given as [14] 

n(M,m,T) = (m ° 4 ~ m)2 +U(*,*,T)-2N f N c J^E p 

_2N f T J ^3 [In (1 + N^e-^ + N c ^e~ 2 ^ + e" 3 ^) 

+ In (1 + N c $e- fiE > + N c ^e~ 2 ^ + e" 3 ^)] . (46) 

Minimizing this thermodynamical potential with respect to $ and $, we obtain equations 

on dtt 

a¥ = °- M = °- (47) 

In the absence of chemical potential, these two equations are identical, and so $ = $ [27]. 
In the same way, Minimizing the thermodynamical potential in Eq. (46) with respect to the 
value of the constituent quark mass m, the gap equation in Eq. (12) can also be obtained. 

In the following, we present our numerical results. First of all, we give our calculated 
values for several characteristic temperatures. These characteristic temperatures include 
the pseudo-transition temperature for chiral crossover, T x , corresponding to the maximum 
of — dm/dT [14,27,29], the pseudo-transition temperature for deconfinement crossover, Tp, 
corresponding to the maximum of d$/dT, the Mott temperature T M for 7r meson, defined 
by 

m n {T M ) = 2m(T M ), (48) 

meaning that the pion can dissociate into a constituent quark and an antiquark above the 
Mott temperature, and the dissociation temperature for a meson TJ" [64], defined by 

m a (T° d )=2m^T° d ). (49) 
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These characteristic temperatures except for T P , can also serve in the conventional NJL 
model. Numerical results for these characteristic temperatures in the PNJL model with two 
Polyakov-loop effective potentials are shown in Table 4. Here, for comparison we also list 
the results in the conventional NJL model. 

Table 4: Several critical temperatures in the conventional NJL model and the PNJL model 
with two Polyakov-loop effective potentials (T = 270 MeV is chosen for these two effective 
potentials). 





T x (MeV) 


T P (MeV) 


T M (MeV) 


T° d (MeV) 


PNJL (Wpoi) 


253.2 


245.4 


264.6 


253.0 


PNJL (Z4n P ) 


245.0 


232.0 


259.6 


246.3 


NJL 


184.4 




201.2 


181.9 



In Fig. 2 we illustrate our calculated results of the masses of it and a mesons, the mass 
of constituent quark, and the Polyakov-loop as functions of the temperature. Fig. 2 shows 
evidently that at a temperature not very high, the masses of the constituent quark, the pion 
and the a mesons maintain the same as the corresponding one at zero temperature. As 
the temperature is around the critical one, these masses vary abruptly. And further, if the 
temperature is very high, the masses of ir and a mesons become degenerate, which indicates 
that the chiral symmetry is restored at high temperature. Such a feature is consistent with 
that given in the framework of Bethe-Salpeter equation combining with the Dyson-Schwinger 
equations (see for example Ref. [65]). We also find that two different Polyakov-loop effective 
potentials do not result in qualitative differences but only slightly quantitative deviations 
as the left panel of Fig. 2 shows. Furthermore, looking through the right panel of Fig. 2, 
we can notice that the chiral phase transition occurs at relatively lower temperature in the 
conventional NJL model. 

In order to compare the obtained results in the PNJL model with those in the conventional 
NJL model more conveniently, we scale the temperature in unit of Mott temperature T M 
and re-display the results in Fig. 3. One can recognize that, in the PNJL model, only 
when the temperature is very near the phase transition temperature, masses of mesons and 
constituent quark begin to deviate from their values at zero temperature obviously. While 
in the conventional NJL model, these masses begin to deviate from their zero-temperature 
values at much lower temperature, about 0.4 T M . This phenomenon can be attributed to the 
fact that, in the low temperature, chiral symmetry is broken and the quark and antiquark are 
in the confined hadronic phase in the PNJL model, contributions from thermal excitations 
of one and two quarks or antiquarks are suppressed as the distribution functions in Eq. (10) 
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Figure 2: Left panel: calculated masses of tt meson, a meson, constituent quark and the 
Polyakov-loop as functions of the temperature. Here, thick curves and thin curves correspond 
to the results with the polynomial Polyakov-loop effective potential U po \, the improved ef- 
fective potential Ui mp , respectively. Right panel: masses of n, a, and m as functions of the 
temperature in the PNJL model with U imp (thick curves) and in the conventional NJL model 
(thin curves). 

and Eq. (11) show when the Polyakov-loop $ approaches zero. This is a manifestation of 
color confinement on the level of statistics. While in the conventional NJL model, due to 
the lack of the appearance of color confinement, contributions from one and two quarks or 
antiquarks become significant even at low temperature, which results in the phenomenon 
mentioned above. 

In Fig. 4 we show the square of the pion-quark coupling strength g% and the pion- 
decay constant f n as functions of the temperature in unit of Mott temperature in the PNJL 
model with polynomial and improved effective potentials and in the conventional NJL model. 
As Eq. (31) shows, when temperature approaches the Mott temperature Tm from below, 
i.e. when the mass of n meson is about twice mass of the constituent quark, iK(m n ) — > 
oo. Therefore, the pion-quark coupling strength and pion-decay constant vanish at Tm, as 
Eq. (28) and Eq. (39) show. Furthermore, One can also find in Fig. 4 that, in the PNJL 
model, g\ and f n almost keep invariant with the increase of the temperature when the 
temperature is not high and these two quantities decrease rapidly only when the temperature 
is above 0.8 T M . While in the conventional NJL model these two quantities begin to decrease 
at about 0.4 Tm- This behavior is also due to the lack of the color confinement in the 
conventional NJL model as the same as the behavior of masses of mesons and constituent 
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Figure 3: Calculated masses of it meson, a meson, and constituent quark as functions of the 
temperature in unit of Mott temperature T M in the PNJL model with U imp (thick lines) and 
in the conventional NJL model (thin lines). 

quark as functions of temperature shown in Fig. 3. 

We have shown above that the Goldberger-Treiman relation and Gell-Mann-Oakes- 
Renner relation at finite temperature are different from those at zero temperature in that a 
factor r defined in Eq. (41) is introduced. The calculated behavior of r as function of temper- 
ature is displayed in Fig. 5 and one can find that, when the temperature is below 0.9 T M , r is 
almost a constant very near 1, which indicates that these two important relations at vacuum 
still serve well in the large region of temperature ~ 0.9 T M . However, when the tempera- 
ture is above 0.9 Tm, r decreases very rapidly and vanishes at T = Tm- In the region of the 
temperature 0.9 T M ~ T M) the Goldberger-Treiman relation and Gell-Mann-Oakes-Renner 
relation at vacuum should be extended to Eq. (40), Eq. (45), respectively. 

4 tt—tt scattering lengths 

The formalism of s-wave n-n scattering lengths at zero temperature in the conventional NJL 
model has been established in Refs. [63, 66, 67], and it has been extended to the case at finite 
temperature by Quack et al. [64]. In this work we follow the notation and calculation given 
in Ref. [63]. The invariant amplitude of n-n scattering has the form: 

( cp c ; dp d \iM \ap a ; bp b ) = iA(s, t, u)5 ab 5 cd + iB(s, t, u)5 ac 5 bd + iC(s, t, u)5 ad 5 bc , (50) 

where a, b, c, and d are the isospin labels, and s, t and u are the Mandelstam variables, 
s — (Pa + Pb) 2 , t — (p a — Pc) 2 and u = [p a — p d ) 2 - The amplitude of definite total isospin /, 
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Figure 4: Left panel: calculated results of the square of the pion-quark coupling strength 
g\ (in Eq. (28)) as a function of the temperature in unit of Mott temperature T M in the 
PNJL and the conventional NJL model. Right panel: calculated results of the pion-decay 
constant /„. (in Eq. (38)) as a function of the temperature in unit of Mott temperature. 
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Figure 5: Calculated factor r defined in Eq. (41) as a function of the temperature in unit of 
Mott temperature in the PNJL and the conventional NJL model. 
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defined by Aj, can be projected out, given by Ref. [63] 

A = 3A + B + C, A 1 = B - C, and A 2 = B + C. (51) 
When the scattering is at the kinematic threshold, we obtain the scattering lengths, i.e. 

ai = Ms = Ami, t = 0, u = 0). (52) 

OZTT 

For simplicity, the pion momenta can be chosen as 

Pa = Pb = Pc = Pd = P, and p 2 = ml, (53) 

which can be verified to fulfill the threshold condition in Eq. (52). To lowest order in 1/N C , 
there are two types of Feynman diagrams contributing to the s-wave tt-tt scattering [63, 66], 
i.e. the box diagram and the a-propagation diagram. Here we also present them in Fig. 6. 
The three diagrams in the first row of Fig. 6 are the box diagrams and the ones in the second 
row are the a-propagation diagrams. 



7T l+p VT 7T 7T 7T l+p VT 




d e f 

Figure 6: Feynman diagrams contributing to the s-wave tt-tt scattering (see also Ref. [63, 66]). 
Here the external momenta for pions are chosen to be the special case in Eq. (53). 

Following the calculation of Ref. [63], we obtain respective amplitude for each diagram 
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in Fig. 6 as 

iM a 



(SabScd + 5 ac 5 bd - 5 ad 5 bc )(-AN c N f g 4 wqq )[I(0) + J(p) - p 2 K(p)] 

= (Sabred + $achd — 5 a dSbc)iT a , 

iM b = (5 ab 5 cd -5 ac 5 bd + 5 ad 5 bc )(-AN c N f gi qq )[I(0) + I(p)-p 2 K(p)} 



(Sab^d - 8 ac 5 bd + 5 ad 5 bc )iT b , 

(SabScd + SacSbd + 5 ad 6 bc )(-8N c N f g* gq ) 

(Sabred + Sacfibd + S ad S bc )iT c , 



1(0) + V -L(p) - 2p 2 K(p) 



(54) 
(55) 

(56) 



5 ab 5 cd (8N c N f gt qq ) 



= 8ab$cdiT d , 



I\p) 



(l - £) 7(2p) + ^I(m w ) 



SMSN c N f gt qq ) m - p : K ^ ]2 



= 8 ac 5 bd iT e , 
iM f = 5 ad 5 bc (8N c N f g 



nqqj 



\m-p^K{p)f 

i(o) + g/w 



= 5 ad S bc iT f . 



(57) 



(58) 



(59) 



From above equations, one can notice T b = T a , and Tf = T e . Substituting Eqs. (54) — (59) 
into Eq. (50), one obtains 



A = 2T a - T c + T d , 
Therefore, employing Eq. (51) we have 



B = C = T c + T e . 



(60) 



A = 6T a -T c + 3T d + 2T e , 

Ai = 0, 

A 2 = 2(T C + T e ). 



(61) 



In Fig. 7 we present our calculated results of the scattering amplitudes T a , T c , T d , and T e 
as functions of the temperature in unit of T M in the PNJL model with polynomial Polyakov- 
loop effective potential and improved effective potential and in the conventional NJL model. 
The results of the conventional NJL model in our work are roughly consistent with those 
given in Ref. [64]. However, there exists a difference which reads that our present calculation 
indicates that the scattering amplitude T a approaches zero at the Mott temperature Tm, the 
calculation in Ref. [64] gives that T a is divergent at T M - Recalling the analysis above, we 
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Figure 7: Calculated scattering amplitudes T a , T c , T d , and T e as functions of the temperature 
in unit of Mott temperature in the PNJL and the conventional NJL models. 

would emphasized that, when the temperature approaches to Tm, iK(m n ) in Eq. (31) and 
iL(m n ) in Eq. (32) are divergent and the degree of divergence of iL{m n ) is higher than that 
of iK(m n ). Substituting the expression of the coupling between n meson and quarks gl qq 
in Eq. (28) into Eq. (54), we find T a oc l/(— iK{m^)) when the temperature is near the 
Mott temperature, therefore, T a approaches zero at T M . In the same way, we find that T d 
approaches zero, T c is divergent, and T e approaches a finite value at the Mott temperature. 
Furthermore, when the temperature is equal to the dissociation temperature of a meson, i.e. 
T = T%, we have m a = 2m n , which results in that the a propagator in diagram d of Fig. 6 
and also the amplitude Td (see Fig. 7) are divergent at TJ. Comparing the results of the 
PNJL model with those of the conventional NJL model, one can also recognize the similar 
behavior as obtained previously, which reads that the T-matrix amplitudes calculated in the 
PNJL model deviate from their values at zero temperature only when the temperature is 
near the critical temperature, while the deviation occurs much earlier in the conventional 
NJL model. Taking amplitude T d for example, since the mass of a meson decreases with 
the increase of the temperature much earlier in the conventional NJL model than that in 
the PNJL model, as Fig. 3 shows, we expect that the divergence in T d also occurs earlier 
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in the conventional NJL model, which is verified in Fig. 7. In addition, we find that the 
value of T%/T M calculated in the conventional NJL model is about 0.90, smaller than 0.96 
in the PNJL model with polynomial effective potential, and 0.95 in the PNJL model with 
improved effective potential. 
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Figure 8: Calculated s-wave ir-ir scattering lengths ao and a 2 as functions of the temperature 
in unit of Mott temperature in the PNJL and the conventional NJL models. 

In Fig. 8, we show the calculated results of the s-wave n-n scattering lengths a and a 2 
as functions of the temperature in unit of Tm in the PNJL model with two Polyakov-loop 
effective potentials and in the conventional NJL model. Since ao and a 2 contain contribution 
from the T- matrix amplitude T c as Eq. (61) shows, they are both divergent at T = T M - 
Furthermore, a also contains T d , so it diverges at T = TJ" as well. At zero temperature, we 
have a = 0.173 and a 2 = —0.045, which are consistent with the Weinberg values (a ) w = 
7ml/(32nf%) = 0.158 and (a 2 ) w = -2ml/(32nf%) = -0.045 [64,66]. On the experimental 
side, the Geneva-Saclay collaboration provided the often quoted values a = 0.26 ± 0.05 and 
a 2 = —0.028 ± 0.012 [68,69], and recent years experiment E865 at Brookhaven National 
Laboratory, USA, has given new values a = 0.203 ± 0.033 ± 0.004 and a 2 = —0.055 ± 
0.023 ± 0.003 [70], and also a = 0.216 ± 0.013 ± 0.004 ± 0.005 [71]. The scattering lengths 
ao and a 2 are almost independent of the temperature until the temperature is increased to 
0.9Tm in the PNJL model, and after that they vary drastically with temperature. While 
in the conventional NJL model a and a 2 begin to vary with temperature at about 0.6T A f 
and the temperature at which ao diverges due to the a-meson dissociation is also lower in 
the conventional NJL model. From the above analysis, one can recognize that the physical 
meaning of the divergence of the s-wave n-n scattering lengths a and a 2 at the Mott 
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temperature T M for ix meson is clear and consistent with that shown in Ref. [64]. At T = T M , 
n meson can dissociate into a constituent quark and an antiquark, the n meson is then 
unbound and its radius becomes infinite. Mathematically, the geometrical size of the pion 
meson can be described by its charge radius r which is related to f w through / 2 oc l/(r 2 ) [64]. 
Employing the Weinberg relations cited above, we have the relation between the scattering 
lengths and the charge radius of n meson as \a\ oc (r 2 ), which clearly indicates that the 
divergence of n-n scattering lengths at the pion Mott temperature is closely related with the 
melting of the pion meson. The relation between the divergence of the s-wave n-n scattering 
lengths and the delocalization of the pion meson at Tm is then confirmed not only in the 
conventional NJL model but also in the PNJL model. As for the divergence of the s-wave 
7r-7r scattering length a at the dissociation temperature for a meson TJ", we should note 
that this divergence corresponds to the situation that the propagator for a meson displayed 
in the part d of Fig. 6 is on shell, which means that the ti-ti scattering in the s channel 
couples resonantly with the a meson field. The divergence (from positive infinite to negative 
infinite) of the a and the mass relation m a = 2m 7r suggest that, in general point of view, a 
very loosely bound state may appear at the dissociation temperature T%. However, detailed 
investigation is required to clarify its mechanism. 



5 Summary 

In summary, we have studied the mesonic excitations at finite temperature in the two flavor 
PNJL model. The masses of n meson and a meson, pion-decay constant, and the pion-quark 
coupling strength at finite temperature are calculated in the PNJL model with two forms 
of Polyakov-loop effective potential. Their variation behaviors with temperature, especially 
when the temperature takes a value near the critical one, are investigated in details. We 
find that the results calculated in the PNJL model are almost independent of the choice of 
the Polyakov-loop effective potential. We also compare our calculated results in the PNJL 
model with those in the conventional NJL model. We find that, since in the PNJL model, 
the Polyakov-loop which is coupled with quarks suppresses the unwanted degrees of freedom 
below the critical temperature, all quantities describing the properties of mesons deviate 
from their values at zero temperature only when the temperature is very near the critical 
temperature, for example the Mott temperature, and they vary with temperature rapidly 
in a very narrow region near the critical one. While in the conventional NJL model, these 
quantities begin to vary with temperature much earlier. Therefore, we conclude that the 
effect of color confinement screens the effect of temperature below the critical temperature. 
Furthermore, we have investigated the Goldberger-Treiman relation and the Gell-Mann- 
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Oakes-Renner relation at finite temperature in the PNJL model, and we find that when 
the temperature is below about 0.9T M , where T M is the Mott temperature for ir meson, 
these two important relations are hardly changed by the effect of temperature. However, the 
Goldberger-Treiman relation and the Gell-Mann-Oakes-Renner relation should be corrected 
once the temperature is in the region of 0.9T M ~ T M - 

In this work, we have also investigated the s-wave n-n scattering lengths at finite tem- 
perature in the PNJL model. The obtained results in the PNJL model are also compared 
with those in the conventional NJL model. We find that scattering length a is divergent 
at Mott temperature for 7r meson, Tm-> and at dissociation temperature for o meson, 
and scattering length a 2 is divergent at T M , which are consistent with the results in the 
conventional NJL model. Due to the effect of color confinement, the dissociation tempera- 
ture for a meson T% calculated in unit of T M in the PNJL model is relatively larger than 
that given in the conventional NJL model. In the same way, the influence of the temper- 
ature on the scattering lengths a and a 2 below the critical temperature is suppressed by 
the color confinement in the PNJL model. In addition, the characteristic of the scattering 
amplitude T a at the Mott temperature calculated in the PNJL model is different from that 
given previously in the conventional NJL model. 
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